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Abstract We analyze the conditions, which guarantee the existence of periodic and 
soliton-like traveling wave solutions in the non-local hydrodynamic model of structured 
media. 

1 Introduction 

In natural science there exists a number of examples of the creation and stable 
evolution of so called coherent states, or spatio-temporal patterns [HOE]. 
The location and specification of such pattens within the confines of adequate 
mathematical models is quite difficult, because they are usually described by 
non-linear partial differential equations. As exceptions, the bi-Hamiltonian 
equations should be mentioned (see e.g. [I], ch. VII), which can be integrated 
completely by applying the inverse-scattering method or other associated 
procedures [5|. Unfortunately, the number of equations which fall under this 
method is scarce. And, what is possibly more important, coherent structures 
occur very often in open dissipative systems which cannot be Hamiltonian in 
principle. In such cases, only a few approaches, the symmetry-based methods 
[61 H] among them, could serve as alternatives to numerical modeling. It is 
worth noticing that in recent years additionally, and seemingly independently 
to the group methods, the so called anzatz-based approach is widely used 
[3 El El EE]. It allows for solutions with pre-set properties (e.g. soliton-like, 
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kink-like or periodic stationary structures or traveling waves) to be obtained, 
but, according to our belief, they are less universal than approaches based on 
self-similarity methods, supported by qualitative analysis. This is because 
in the latter case one can analyze not only particular self-similar solutions, 
which fortunately can be expressed analytically, but the whole family with 
the given symmetry. 

Let's now focus our attention on a very important factor which is inten- 
sively used in the article. Dissipative models of real processes very often 
turn out to be close (to some extent) to Hamiltonian models. Such a situa- 
tion occurs for example in cases when dissipation, relaxation, heterogeneity 
and similar effects which ruin the regularity of the model, display a small- 
parameter description. In such while searching for dissipative struc- 
tures, not only qualitative analysis can be applied, but also one can use the 
whole power of Hamilton's method. Apart from this, to research a system 
with small parameters, the ideas and methods of the concept of approximate 
symmetry can be used, and physically justified small modules, regardless of 
whether or not they break the symmetry of an unperturbed model, may be 
taken into consideration. 

This article is concerned with the search for wave structures (particularly 
soliton-like and periodic regimes) in a hydrodynamic model, taking into ac- 
count the effects of spatio-temporal non-localities. Using the afore-mentioned 
methods, the criteria for the existence of the wave patterns will be formulated 
in a wide set of parameter values in those cases where the model under con- 
sideration is in some way equivalent to a Hamiltonian one, and in the case, 
when dissipative effects (associated with relaxation) and the mass forces are 
taken into account. 

The structure of the article is as follows. In section 2 a hydro dynamic-type 
model of a non-local medium is formulated, and a group-theory reduction to 
the system of ODEs, describing a family of approximately invariant travel- 
ing wave solutions is performed. In section 3 the existence of one-parametric 
family of periodic and soliton-like solutions is shown when the temporal re- 
laxation and mass forces are ignored. In section 4, based on a generalization 
of Melnikov's method, and the essential use of the previously obtained exact 
solution, the existence of soliton-like regimes is demonstrated when small 
mass force and relaxation effects are incorporated. 
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2 Nonlocal hydrodynamic-type system and 
its symmetry reduction 

Below we introduce a modelling system, taking into account the non-local 
effects. These effects are manifested e.g. when an impulse loading is applied 
to media, possessing an internal structure on mesoscale [HI H21 E3] • Descrip- 
tion of the non-linear waves propagation in such media depends in essential 
way on the ratio of a characteristic size d of elements of the medium struc- 
ture to a characteristic length A of the wave pack. If d /X is 0(1), then the 
basic concepts of continual media are completely false and one should use 
the description, based, e.g. on the element dynamics methods [H] . The clas- 
sical continuum mechanics equations applications is justified rather in those 
cases, when d /X <1, and the discreteness of the matter could be completely 
ignored. 

The models studied in this work apply when the ratio d/X is much less 
that unity and therefore the continual approach is still valid, but it is not 
as small that we can ignore the presence of the internal structure. As it 
have been shown in a number of papers (see e.g. [15]), in the long wave 
approximation the balance equations for mass and momentum retain their 
classical form, which in the one-dimensional case can be written as follows: 



u t + p x = %/p, ,* 
Pt + P 2 u x = 0. 1 ' 

Here u denotes mass velocity, p is pressure, p is density, t is time, x is 
mass (lagrangean) coordinate, is mass force, lower indices denote partial 
derivatives with respect to subsequent variables. Thus, the whole informa- 
tion about the presence of structure in this approximation is contained in a 
dynamic equation of state (DES), which should be incorporated to system 
(JIJ in order to make it closed. 

Generally speaking, DES for multi-component structured media, mani- 
festing the non-local features, takes the form of integral equation [THJ ITT] , 
linking generalized thermodynamical flow J and generalized thermodynami- 
cal force X, causing this flow: 



J 



dt'. (2) 



K(t, t'; x, x')X(t', x') dx' 

Here K(t,t'; kernel, taking into account nonlocal effects. Function 

K can be calculated, in principle, by solving dynamic problem of struc- 
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ture's elements interaction, however such calculations are extremely difficult. 
Therefore in practice one uses, as a role, some model kernel, describing well 
enough the main properties of the non-local effects and, in particular, the 
fact that these effects vanish rapidly as \t — 1'\ and \x — x'\ grow. This prop- 
erty could be used in order to pass from the integral equation (T5]) to a pure 
differential equation. 

One of the simplest state equations accounting for the effects of spatial 
nonlocality takes the form: 



p = a 



K x {t, t') 



+00 



K 2 (x, x') p n (t',x')dx' 



dt\ 



(3) 



where K 2 (x, x') = exp [— (x — x') 2 /l 2 } (cf. with [13J). Using the fact that 
the function exp [— {x — x') 2 /l 2 } extremely quickly approaches zero as \x — x'\ 
grows, we substitute the function p n (t',x') by the first three terms of its 
decomposition into the power series: 



p n (t\x') = p n (t\x) + [p n (t\x)] x 



X 



0{\x-x'\ 3 ), 



1! ■ L ' v ' " 2! 
obtaining this way the following approximate flow-force relation: 

ft 



p = a Ki(t, t') L (p, p x , p, 

J —00 



, cj dt , 



(4) 



where 



L[p, p x , p xx ] = c p n (t',x) + c 2 [p n (t',x)} 



+00 



Cq — I I e T dr = 



(■2 



2 



+00 



r 2 e- T 



XX ' 

dr 



Now we are going to analyze different function K\ (t, t'), responsible for 
the relaxing effects inside the elements of the internal structure. A medium 
with one relaxing component is usually described by the kernel 



K x (t, t') = -T 1 exp 



t-t' 



T 



(5) 



With such a kernel the flow-force relation (j4j) can be expressed as a first 
order PDE. In fact, inserting the kernel (jHJ) into the equation (j4j) and then 
differentiating it with respect to t, we get the following DES: 



rp t +P = &L[p, p x , px 



(6) 
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Assuming that 



K\ (t, t') = at exp 



t-t' 



+ a 2 exp 



t-t' 

T 2 



we can obtain the DES for structured media with two relaxing components. 
In order to obtain a pure differential flow-force relation, we should use p, p t 
and pu- Taking their linear combination with properly chosen coefficients, we 
obtain the following DES: 



hptt + rp t +p = aL[p, p x , p xx ] + pL t [p, p x , p xx ) 



(7) 



where 



h = ri r 2 , t = t 1 + t 2 , a = a [a x t x + a 2 r 2 ) , p = a (a x + a 2 ) n r 2 . 

It is worth noting that a DES formally identical with ((7j) can be obtained 
from the equation ((3]), using in the temporal part of the kernel function 



K\ (t, t') = exp 



t-t' 

Tl 



COS 



t-t' 



r 2 



+ Vo 



which contains, beside the exponentially decaying, the oscillating term. In 
fact, using the similar combination as in the previous case, we get the equa- 
tion, formally identical with the Eq. ([7]): 



hptt + Tp t +p = aL[p, p x , p xx ] + pL t [p, p x , p a 



(8) 



with 



h 



n 



7j T = 7, a = h cosv9 7, H 



7 

2r 2 



(n cos ip - t 2 sin (p ) , 



where 7 = 2t x rf/ (rf + r|) . In order to maintain the physical meaning of 
pressure and density we introduce the variables p = p—p , p = p — p , where 
p > and p > are some constant equilibrium values of the parameters. 

The difference between ([7]) and ([8]) arises from the fact that the coeffi- 
cients of these formally coinciding equations belong to distinct domains of 
the parameter space. For example, the coefficients h and r from ([7]) satisfy 
the relation 

Ah < t 2 
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provided that T\ ^ r 2 while for the h and f the opposite inequality holds. 
As it is shown in [TBI [T9~] . the properties of traveling wave solutions to a 
non-local hydro dynamic-type model depend in essential way on the values of 
the parameters. 

It is worth noting that DES very similar to © have been obtained in [201 
|2"T] on the basis of the phenomenological non-equilibrium thermodynamics 
formalism. We prefer to present the approach based on the formula ([2]) and 
the modeling kernels of non-localities, since it is less cumbersome than that 
using the thermodynamical methods. 

In this work we concentrate on the study of the system of balance equa- 
tions flU, closed by the DES (JgJ: 

u t +p x = Ss/p, 

Pt + P 2 u x = 0, (9) 
rpt +P = £~ 2 P U+2 + o[p v+l p xx + {y + l)(f (p x )% 

where v = n — 2, = cq a {y + 2) f , a = c 2 o {y + 2) f. 

In what follows we put in <$f) f = er, Q/p = ef(p) with |e| 1. In- 
troduction of a small parameter enables us to employ the approximate sym- 
metry technique and avoid difficulties arising from the fact that the external 
force, as well as the relaxing terms, can destroy the scaling symmetry of the 
problem. 

It is easy to verify by the straightforward checking, that for e = the 
system admits the following Lie symmetry group [U E] generators: 

^o = £, A = &, Z = -^ l tf- t+ uf u + ^ lP f p + ^pf p . (10) 

For small e, the system ([9]) admits the operator 

l = P + D?i+4 (11) 

in the sense of approximate symmetry [22J . Let us note that the first equation 
of the system (Q admits the operator Z in generally accepted sense when 
f(p) = ap u+2 . Nevertheless the set ffTU]) is not admitted by the system (JH]), 
because of the time derivative term contained in the DES. 

A passage to the self-similar variables enabling to factorize the system 
Qj, is based on the operator (TTTj) . The characteristic system corresponding 
to this operator, up to 0(e 2 ), can be presented as follows: 

(v + 3 \ , dx du dp dp 
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Solving the system (j!2p . and expressing the initial variables in terms of its 
first integrals, one can construct the following ansatz: 

u=(l + e£t)tf(n), p= (l + ^e^)n(fi), p= (i + ^e^t)^), 

(13) 

where 

is new independent variable. 

The approximate self-similar reduction is performed in several steps. In- 
serting ansatz ffl3|) into the second equation of the system (j5J), we obtain a 
first order differential equation admitting the separation of variables: 

R 2 if - D R + — — e £ R = 0(e 2 ). 

Integration of this equation gives us the approximate first integral 

u ^= c ^-^Sm +0(t2) - (14) 

In what follows we assume that C\ = D/Rx, where < R\ — const. For 
e = such a choice immediately leads to the asymptotics 

lim u(t, x) = 0, lim p(t, x) = R\. (15) 

x— >+oo x— >+oo 

Inserting ansatz (|T5|) into the first equation of the system (jSj), and using 
Eq. ([Hj), we obtain the equation 

n - D (§R - I + 1 £ |) + e( ( Cl - g - « S(R) - 0(f). (16) 
Now we introduce new function: 

D 2 

G = n + — . (17) 

Taking derivative of (JT7J) with respect to Q and employing (|16p . we obtain 
the equation 



G 



0(e 2 ). (18) 
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Inserting ansats (TT5]) into the third equation, we get a second-order ODE. 
Excluding II from this equation by means of the formulae f)16p . (IT7|) . and 
introducing new variable R = Y we finally obtain a closed system which, up 
to 0(e 2 ), takes the form 

R = Y, 



a R v+2 Y = G R- 



D 2 + 



ft 



R v+3 + a{v + l)R v+l Y 



u + 2 



D 3 

-re—Y, (19) 



G 



vD 



1 + vR, 

In the following sections we analyze in detail the system (TT9~j) for e = and 
formulate the conditions of the existence of periodic solutions as well as the 
homoclinic regimes, corresponding to soliton-like travelling wave solutions of 
the initial system of PDEs. 



3 Qualitative analysis of the reduced system 
in the case when e = 

Assuming that e = 0, we immediately get that G = Gq = const, and the 
system (TT9~j) reduces to 



dR _ y 

H = (aR^ 2 )- 1 {GR-[D 2 + £- 2 R u+3 + <t(u + l)R v+1 Y 2 ] } , (20) 



where u = x — Dt. Incorporating the conditions ([To]) , we express Go as 
follows: 

G ^C + ^ +2 ' (21) 

Dividing the second equation of the system fl20|) by the first one, and in- 
troducing new variable Z = Y 2 = (dR/du) 2 , we get, after some algebraic 
manipulation, the linear inhomogeneous equation 

Z\R) + 2[(u + 1) R^ZiR) = 2[GR-D 2 - f3R u+3 /(v + 2)] /(aR u+2 ). 

(22) 
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Solving this equation with respect to Z = Z(R) and next integrating the 
equation obtained after the substitution Z = (dR/du) 2 , we get the following 
quadrature 

f ±^R^dR 

u - u = / (23) 



(2+z/) (l+i/) > J (2+i/)2 

The direct analysis of the formula (12"3"|) is rather difficult. To realize what 
sort of solutions we deal with, the methods of qualitative analysis can be 
applied to the dynamical system (|20|) . It is evident, that all isolated critical 
points of the system (1201) are located on the horizontal axis OR. They are 
determined by solutions of the algebraic equation 

P(R) = -J^R U+3 ~ GR + D 2 = 0. (24) 

As can be easily seen, one of the roots of equation ( 124"]) coincides with R\. 
Location of the second real root depends on the relations between the pa- 
rameters. If v + 3 > 1, and D 2 satisfies the inequality 

D 2 > D 2 cr = PR v +\ (25) 

then there exists the second critical point R2 > R\- Moreover, if < v is a 
natural number, then the polynomial P(R) has the representation 

P(R) = (R-R 1 )(R-R 2 )^(R), (26) 

where 

^, R . = p[r»+ 1 +r»{R 2-Ri)+---+r{r» 2 -r\)HR2 +1 -K +1 )\ 



{v+2)(R 2 -R 1 ) 

Note that ^(R) is positive, for R > 0. It is a direct consequence of the 
existence of representation (126!) valid for any natural v. But this is also true 
for any v > —2, or, in other words, whenever the function R u+3 is concave 
for positive R. 

Analysis of system's ( )20l) linearization matrix 



M(R U 0) 



1 

(aR^Y^iR^R, -Ri) 



1,2, j^i 

(27) 
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shows, that the critical points A 1 (i? 1 ,0) is a saddle, while the critical point 
A 2 (i?2,0) is a center. Thus, the system (l2"Uj) has only such critical points, 
which are characteristic to a Hamiltonian system. This circumstance suggests 
that there could exist a Hamiltonian system equivalent to the system (120]) . 
The Hamiltonian function would help us to make a complete study of the 
phase portrait of system f[2T)j) . and find out homoclinic trajectories, which 
correspond to a soliton-like wave packs, providing that such trajectories do 
exist. 

If we introduce a new independent variable T, obeying the equation = 
aR u+2 (p(R,Y)-£^, then we can write down the system (1201) as follows: 

§ =aR^Y<p(R,Y), (9R , 
% = <j>{R, Y) (GR - [D 2 + £- 2 R v+z + o(y + l)R» +l Y 2 ] ) . ^ 

Here <f)(R, Y) is a function, that is to be chosen in such a way that the system 
(I28p be Hamiltonian. So we are looking for a function H(R,Y), satisfying 
the system 

^ = aR^ 2 4>(R,Y), 

H = -<f>{R, Y) [GR - \d 2 + -^—R u+Z + o{y + l)^ +1 y 2 l ) . 
oR y z/ + 2 J J 

Equating mixed derivatives of H, we obtain the characteristic system 
dR dY dd> 



(29) 



aR"+ 2 Y GR-[D 2 + £-R»+ 3 + a(v + l)R»+ l Y 2 ] vaYRy+\ 
The general solution of system ( 129]) can be presented in the form 

where $(•) is an arbitrary function of the variable 

p = R"+i {F 2 $ + D 2 a- 1 R-( V+ V log R + R 2 (5/[a(u + 2){v + 3)] - Ga^R^} 
Putting = 2R U , we can easily restore the Hamiltonian function: 

H = 2D 2 — + $ p2(„+2) + aY 2 R 2 ^ _ 2G ^1 (30) 
+ l (z/ + 2) 2 z/ + 2 v ; 

As is well-known, the function H is constant on the phase trajectories of 
the both system f[2"0"j) and (1251) . and since the integrating multiplier = 2i?^, 

10 



appearing in the formula (1281) is positive for R > 0, then the phase portraits 
of the systems (I2UI) and (I28p geometrically are identical in the right part of 
the phase plane (R, Y). Thus, all the statements concerning the geometry 
of the phase trajectories of the system fl28|) lying in the right half-plane is is 
valid to the corresponding solutions of the system (120]) . 

Using the linear analysis, we have already shown that the critical point 
A 2 (R2, 0) is a center. Stated above relations between systems f l2"Ul) and fl2"51) 
enable to conclude that this point does not change when the nonlinear terms 
are added. This means that the critical point A 2 (_R 2 , 0) is surrounded by the 
closed trajectories, and, hence, the unperturbed source system (J9]) possesses 
a one-parameter family of periodic solutions. If the right branches of the 
separatrices of the saddle point A\{Ri, 0) go to infinity (the stable branch 
W s when t — > — oo and the unstable branch W u when t — > +oo), then 
the domain of finite periodic motions is unlimited. Another possibility is 
connected with the existence of the limiting trajectory, bi-asymptotic to the 
saddle. In this case the domain of periodic solutions is bounded, and the 
source system, in addition to a one-parametric family of periodic solutions, 
possesses localized soliton-like regimes, corresponding to the homoclinic loop. 
To answer the question on which of the above mentioned possibilities takes 
place, the behavior of the saddle separatices, lying to the right from the line 
R = Ri should be analyzed. We obtain the equation for saddle separatices 
by putting H = H(R 1 ,0) = Hi in the left hand side of the equation (IBTJj) . 
and solving it with respect to Y: 



It is evident from equation (13T1) , that incoming and outgoing separatrices are 
symmetrical with respect to OR axis. Therefore we can restrict our analysis, 
e.g., to the upper separatrix Y + . First of all, let us note, that separatrix Y + 
forms a positive angle with the OR axis in the point (R\, 0): 



The above formula arises from the linear analysis of the system ( 1201) in critical 
point Ai(Ri, 0). So Y + (R) is increasing when R — R\ is small and positive. 
On the other hand, the function 




(31) 




Q(R) = #i + 2G 



R u+2 



- 2D 



Z/ + 2 



is + 1 (z/ + 2) 2 
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appearing in the RHS of the formula (IHTl) . tends to — oo as R — > +00, because 
the coefficient at the highest order monomial R 2 ( u+2 ^ is negative, while the 
index z/ + 2 is assumed to be positive. Therefore the function Q(R) intersects 
the open set R > R\ of the OR axis at least once. Let us denote the point 
of the first intersection by R%, and let us assume that R3 > R 2 . If, with 
this assumption, we were able to prove that Y± (R) form the right angle with 
the OR axis at the point R%, then we would have the evidence of tangent 
intersection of the stable and unstable saddle separatrices. 

We begin with the note that lim^^- Q(R) = +0. Calculating derivative 
of Q(R) we get: 

Q'(R) = -2R U (^-J^R U+3 ~GR + D 2 ^j = -2R U P(R). (32) 

It follows from the decomposition ( 126]) . that Q'(R) < when R > R 2 . 
Therefore 

dY RQ'(R) — 2{y + l)Q(R) 

hm — — = ■ = —00. 

Rr+Ri dR 2^QR^+ 2 

So, to complete the prof, we must show that the inequality R 3 > R 2 is true. 
Supposing that the inequalities Ri < R 3 < R 2 take place, we obtain from 
the equations f )32|) . f l26|) that lim^^- Y+(R) = +00. On the other hand, 

the function Y + (R) approaches zero remaining positive as R — > R% . But 
such behavior is impossible for any function, being regular inside the interval 
(Ri, R3). The case R3 = R 2 should also be excluded, because the critical 
point A 2 (R 2 , 0) is a center. The result obtained can be formulated as follows. 

Theorem. Ifu > —2 and D 2 > (3 R± +3 , then the system [2D\I possesses a 
one parameter family of periodic solutions, localized around the critical point 
A 2 (R 2 , 0) in a bounded set M. The boundary of this set is formed by the 
homoclinic intersection of separatrices of the saddle point A\ {R\, 0). 

Thus the unperturbed source system (EJ possesses periodic and soliton- 
like invariant solutions. Let us note in conclusion that for some special case 
the integral at the RHS of the formula (|23|) or, what is the same, at the RHS 
of the formula 

T ~ T ° = ± I 2^R^jH 1 + 2G§^-2D^-0^ <33) 
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can be calculated explicitly: 



T 



"7\/2 



7 In 



R-l 



.3- R + V7 -2R- R 2 



+2 V^ln 



7 - R + V7V7 -2R- R 2 



R 



(34) 



This solution corresponds to the following values of the parameters: D 
1 = R 1 = a = 1, = 1/2, i/ = 0, G = 5/4, and T = 0. 



4 Homoclinic solutions of perturbed system 



The purpose of this section is to analyze whether and when the perturbed 
system ffT9l) possesses the homoclinic solution. Our analysis is based on the 
generalized Melnikov theory presented in [23]. The theory is based in essential 
way on the knowledge of the homoclinic solution of unperturbed system, 
which enables to measure the distance between the stable and unstable saddle 
separatrices for small e. For this reason we specify the parameters in such a 
way that they fit the exact solution (|34|) . 

Since we are going to use the formalism developed in |23j, we should pass 
in (|T9|) to the independent variable T, for which the unperturbed system 
becomes Hamiltonian. The standard representation for our system in this 
case will be as follows: 



dX 
dT 



JD X H(X, G) + eF(X, G), 



(35) 



where X = colon (R, Y), D x = colon ^?}, < e< 1, 

F(X, G) = colon(-2D 3 Y/R, 0), 

1 - 



L{X, G) = 2 



D 
R 



while 



J 



In what follows we put f(R) = a R 2 + b. 

So we are going to measure the distance between the saddle separatrices 
of the perturbed system fl35l) . The problem we deal with differs from the clas- 
sical one [2U [22] since in the unperturbed case the Hamiltonian f[3"Ul) depends 
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on an extra variable G, which plays the role of a parameter. We assume that 
G belongs to an open set U, containing G = D 2 / R x + R{ +2 /(v + 2). 
Following [23], we introduce the set 

M = {(X, G) G M 2 x R 1 : X = 7 (G), where j(G) solves 

D x H(-y(G), G) = and det[D 2 x H(^(G), G)] ^ 0,VG e U}. (36) 

We denote by IV s (Ai) and W U (M) the 1 + 1- dimensional stable and unstable 
manifolds of Ai. It is obvious that the manifolds W s (Ai) and W u (Ai) 
coincide as e = (in this case W S (M) = W U (M) = T). For sufficiently 
small e, Ai is transformed into the locally invariant manifold W £ (Ai). Up 
to 0(e 2 ), the dynamics on W £ (Ai) is governed by the equation 



For b = 2£— a , W e (Ai) possesses the equilibrium point (1, 0, 5/4) . Analysis 
of the linear part of equation fl37|) shows, that this point is stable as a > £/2 
and unstable as a < £/2. 

For nonzero e, the invariant manifolds W s (Ai) and W u (Ai) are trans- 
formed into the (locally) invariant manifolds W^(Ai) and W^(Ai), which 
do not coincide. Yet at certain conditions the manifolds can still have the 
points of intersection, different from W e (A4). To state these conditions, a 
distance separating W*(A4) and W^(Ai) is measured in vicinity of a point 
peT, lying in the section G\ = 5/4. Up to 0(e), this distance is equal to 
the Melnikov integral ([23], Ch. IV): 



where X Gl (T) = (R Gl (T),Y Gl (T)) is the homoclinic trajectory of the un- 
perturbed system on the G = G\ level corresponding to the hyperbolic fixed 
point of the vector field on M. £ . The dependence on (X Gl (T), G\, T) means 
in other words that the integration is carried out along the unperturbed 
trajectory V. Formula (134"]) describes this trajectory as an implicit function 
T = T(i? Gl ). 



G = eL{ 1 {G) i G). 



(37) 




+ (D X H,(D G JD X H) L(X,G)dT) (X^(T),G u T)dT, (38) 
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Let us begin the calculation of Melnikov integral. Below we present in 
explicit form the terms from the expression 



D X H = 
JD X H = 
D G JD X H 



'2Y 2 R + 2 + lR 3 -2RG 
2YR 2 



1 
-1 

" ' 

2R 



, (D X H) 



2YR 2 



-2Y 2 R-2-\R z + 2RG 



J LdT = J 2R 2 [a(R Gl (t)) 2 + £ 
Let's replace T with ( j34"j) . Then: 

1 



R°i(t) 



}dT. 



dT 



R 2 



2R 2 \l I - 2R - |i? 4 + \R 2 



dR. 



Finally 



J LdT = J(aR 2 + i-a-^)J-_ 



R 2 



2 R - + \R? 



dR 



V2a [1 - R Gl (T)] y/7-2R G *(T) - (R Gl (T)) 2 + 



2V2{£ + 4a) arcsm 



i? Gl (T) + 1 
2^2 



So the Melnikov integral will be of the form: 

M Gl (a,0 = h-h, 

where 

h = J AR Gl (T) \Y Gl (T)] 2 dT, 

T(G) 

I 2 = 8V2 J R Gl (T) [R Gl (T)] 3 Y Gl (T) { | [l - R Gl (T)] ^6 - [1 + (T)f 



r(G) 



^ . i \ . i? Gl (T) + l 

(£ + 4 a) arcsm = 

2^2 
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where T(Gi) is the unperturbed homoclinic trajectory on the G\ level. 
The computation of 1\ is based on the Green's theorem. As 

Y Gu n _ 1 dR Gl {T) 
K ' 2(R G ^{T)) 2 dT ' 

we have: 

-h= f -A(Y Gl (T)) 2 R Gl (T)dT = I -2^-dR + 0dY 
Jr(G!) /r(Gi) R 

Next, 



rs , (R-l)^/7-2R-R 2 
Y /■2V2-1 f± 2V2R 2 

-2-dR + 0dY = / / -dYdR 



r 



R h I (R-iWr-2R-R 2 R 

2-/2R 



= 1 2(7? 1] ^L R2 2R —dR = !V2(-7+V7hi(8+3V7)) = 0,262789. 

On calculating the J 2 , we use the fact that Y Gl (T) = 2(i? Gl (T))~ 2 dH ^ (T) : 



/ r(Gi) 8y2( J R Gl (T)) 3 r G 1 (T)(-ia( J R G 1 (T) - 1)^1 -2B&{T) - [^(T)] 2 + 
+(£ + 4a) arcsin( fiG ^ )+1 ))ciT = 

= 2 f^ 1 4V2R(-±a(R - 1)V7-2R-R 2 + (£ + 4a) arcsin(|Ji))ciR = 

= -8^-3 + vr)a + (4^ + (-8 + 6y/2)ir)(^ + 4a) = 
= 2 [-1, 60194a + 7, 18141(f + 4a)] . 

(39) 

So, finally 

M Gl (a,0 = 0,262789 -2- [1.60194a + 7, 18141(f + 4a)] , (40) 

amd the Melnikov's integral vanishes if 

f = -0,03659 - 3, 7769a. (41) 

Then, for e sufficiently small, stable and unstable manifold intersect near the 
values of (a, £), defined by equation 

16 



0.3 ; 

o.2; 
o.i 1 



-0.10 



-0.05 



-0.1 : 

-o.2; 
-o.3; 
-o,4; 



0.05 



0.10 



Figure 1: The points (a, £), corresponding to the homoclinic intersection 
(boldface dots) on the background of the straight line (T4"Tj) 



a 


-0.09 


-0.07 


-0.05 


-0.04 


-0.025 


-0.01 


5 


-0.0214 


-0.0186 


-0.01567 


-0.014175 


-0.0115475 


-0.009478 


a 


0.01 


0.025 


0.04 


0.06 


0.08 


0.1 


5 


-0.0065 


-0.0041 


-0.0017 


0.000125 


0.0052 


0.0089 



The results obtained was verified by direct numerical simulation. To 
gain the locus of the points in the parameter plane (a, £) corresponding 
to the homoclinic intersection, we used the following modification of the 
approximate formula (T4"Tj) obtained with the help of the Melnikov method: 



£= -0,03659 - 3, 7769 a + 5. 

The values of the parameters at which the homoclinic intersection is observed 
are presented in the table shown above. The dynamical system (fl9|) was 
numerically integrated under the following values of the parameters: e = 
-0.001, D = l = R 1 = a = l, p = 1/2, v = 0, G = 5/4. The points of 
the parametric plane (a, £) corresponding to the homoclinic intersection are 
shown in figure [U together with the straight line given by the formula (1411) . 
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5 Concluding remarks 



We have considered a non-local hydrodynamic-type model of structured me- 
dia and studied a family of traveling wave (TW) solutions. Our analysis 
shows that in the case of a pure spatial non-locality and the absence of mass 
forces a reduced ODE system is equivalent to a Hamiltonian one. It is worth 
noting that the initial system of PDEs is not Hamiltonian for the physically 
justified values of the parameters. 

Equivalence of the reduced system to a Hamiltonian one is employed in 
this work to obtain rigorous proof of the existence of a one-parameter family 
of periodic and soliton-like TW solutions. At certain values of the parameters 
these solutions possess an analytic representation, while in general they are 
described in terms of Jacobi elliptic functions. 

Having established a soliton-like solution to the rigorously reducible model, 
it is possible to analyze the soliton-like solutions' existence in general case, 
when the exact symmetry of the modeling system has broken down. Us- 
ing, in this case, the ideas of the approximate symmetry and the generalized 
Melnikov method, we've stated in explicit form the conditions of the saddle 
separatrices intersection in cases when the small mass force and the effect 
of temporal non-locality are incorporated. As there arise from the general 
considerations presented in [23], intersections of stable and unstable saddle 
separatrices are not transversal for the given type of perturbation and this is 
backed by the results of our numerical experiments, that enable us to reveal 
the loci of the saddle separatrices intersection with no signs of a "homoclinic 
blow-up" observed. 

In view of this, it is interesting to investigate what type of temporal non- 
locality and external force would lead to the transversal intersection and, as 
a consequence, to drastic changes in the phase portrait of the reduced sys- 
tem. Our preliminary analysis shows that the transversal intersection can 
take place in the case of periodic mass force and the DES (J7|), describing the 
structured media with two relaxing processes. In fact, the second relaxing 
process incorporation causes the dynamic equation of state to include as a 
singular perturbation the second derivative with respect to time. An extra 
additional variable, playing the same role as G, will then appear in the dy- 
namical system obtained via the approximate symmetry reduction. Let us 
note in conclusion, that the role of spatial non-locality on the wave patterns' 
formation and evolution have been discussed in [2T]. In particular it was 
shown in this work that chaotic wave patterns can be gained via homoclinic 
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bifurcation in the hydro dynamic-type model accounting for the effects of spa- 
tial non-localities and where, as in the DES (J7J), the term p tt is incorporated 
small addend. 
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